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It is generally believed that neutrino mass hierarchy can be determined at a long baseline exper- 
iment, often using accelerator neutrino beams. Reactor neutrino experiments at an intermediate 
baseline have the capability to distinguish normal or inverted hierarchy. Recently it has been 
demonstrated that the mass hierarchy could possibly be identified using Fourier transform to the 
L/E spectrum if the mixing angle sin (2^13) > 0.02. In this study a more sensitive Fourier anal- 
ysis is introduced. We found that an ideal detector at an intermediate baseline (~ 60 km) could 
identify the mass hierarchy for a mixing angle sin'^(26'i3) > 0.005, without requirements on accurate 
information of reactor neutrino spectra and the value of Am32. 



PACS numbers: 13.15.-|-g, 14.60.Pq, 14.60.Lm 



Recent results from solar, atmospheric, reactor and ac- 
celerator neutrino experiments all show that neutrinos 
are massive and they can oscillate from one type to an- 
other. Among all the six mixing parameters, three of 
them are known, two unknowns, and one of them, the 
mass-squared difference AtoIj: defined as to| — m|, is 
only known to be lAmj^] = (2.43±0.13) x IQ-^eV^ (68% 
C.L.) from accelerator neutrino experiments ^1]. The 
question, if the mass hierarchy is normal {Am^2 > 0) 
or inverted (Ato§2 < 0), is not known now but is funda- 
mental to particle physics. 

For normal hierarchy (NH) or inverted hierarchy (IH) , 
the neutrino mass-squared difference has the following 
relations: 
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In principle, the mass hierarchy can be determined by 
precision measurements of jArngj^j and |Am32|. In fact it 



is extremely difficult since Am^i is only '^ 3% of lAm^jl 
hence |Am32| and |A7tt,§]^| have to be measured with a 
precision much better than 3 %. 

Effects of mass hierarchy can be amplified by matter 
effects if the baseline is large enough, say several hun- 
dreds to thousands of kilometers. Such experiments often 
need accelerator-based neutrino beams and huge detec- 
tors. Proposals such as T2K 0, Q, Nova [1, [1 [1 and 
T2KK 6] have mass hierarchy sensitivity in the v^ —^ v^ 
channel if 6*13 is large enough (i.e. sin^(26'i3) > 0.03). In 
addition, they are affected by the {Snp, signjArnj^)) de- 
generacy 0,13. At a magic baseline @,[l3l5 L ^ 7000 km, 
the degeneracy can be canceled but it requires a very in- 
tensive source such as a neutrino factory or a beta-beam 
which will not be available in the near future. A method 
using atmospheric neutrinos [111 . Il2| with a baseline of 
L ^ 10'* km and the neutrino energy oi E ^ 1 GeV is 
sensitive to mass hierarchy for very small or even null 
value of 013, if the measurement precision of |A?7i|2| is 
better than 2%. 

Method using reactor neutrino based intermediate 
baseline (40 — 65 km) experiments has been explored 



based on precision measurement of distortions of the en- 
ergy spectrum due to non-zero 6*13 [ij, [l^] . Recently, a 
study 15] shows a new method to distinguish normal or 
inverted hierarchy after a Fourier transform of the L/E 
spectrum of reactor neutrinos. It is observed that the 
Fourier power spectrum has a small shoulder next to the 
main peak, and their relative position can be used to 
determine the mass hierarchy. Afilter method is used 
to improve the sensitivity to the mass hierarchy up to 
sin (26'i3) > 0.02, if A77i§2 is known a priori. Compar- 
ing to a normal L/E analysis, the Fourier analysis natu- 
rally separates the mass hierarchy information from un- 
certainties of the reactor neutrino spectra and other mix- 
ing parameters, which is critical for very small sin (2^13) 
oscillations. 

In this paper, we report that if a proper Fourier trans- 
form is applied and if all information is fully utilized, 
the capability of an intermediate baseline reactor exper- 
iment to determine the neutrino mass hierarchy can be 
improved for a smaller mixing angle ^13 without know- 
ing Arn^2 0- priori. In the following, we will use a reactor 
neutrino spectrum to illustrate the method, but such a 
method can be generalized to other experiments. 

For a reactor neutrino experiment, the observed neu- 
trino spectrum at a baseline L, F{L/E), can be written 
as 

F{L/E)^cp{E)a{E)P,,{L/E) 

where E is the electron antineutrino {pe) energy, (j){E) is 
the flux of Vf, from the reactor, a{E) is the interaction 
cross section of Fe with matter, and Pee{L/E) is the Ve 
survival probability. 

The j/g flux (j){E) from the reactor can be parameter- 
ized as [l6J . 
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where four exponential terms are contributions from iso- 
topes 235U, 23^Pu, 23^U and 24ipu in the reactor fuel, 
respectively. 



parameter 


best fit 


2a 


3cr 


Amii[10-5eV2] 


7.6 


7.3-8.1 


7.1-8.3 


|Ami2|[10-'eV2] 


2.4 


2.1-2.7 


2.0-2.8 


sin^ 6I12 


0.32 


0.28-0.37 


0.26-0.40 


sin^ 623 


0.50 


0.38-0.63 


0.34-0.67 


sin^ 6*13 


0.007 


< 0.033 


< 0.050 



TABLE I: Neutrino mixing parameters from a global fit, up- 
dated in 2007, as the inputs to this study. 



The leading-order expression for the cross section [17 
of inverse-/? decay ( Fe + p — > e^ + n ) is 



a^'^^ = 0.0952 X 10-*^cm^(^^"')pi"'VlMeV^) (3) 

where Ee — E^ — {Mn — Mp) is the positron energy 
when neutron recoil energy is neglected, and pe is the 
positron momentum. The survival probability of Ve can 
be expressed as |18j 
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where Ay — 1.27AmfjL/E, Amf, is the neutrino mass- 
squared difference (mf — rn^) in eV^, % is the neutrino 
mixing angle, L is the baseline from reactor to V^ detector 
in meters, and E is the T7e energy in MeV. 

Pee[L/E) has three oscillation components, P21, P31 
and P32, corresponding to three oscillation frequencies 
in L/E space, which are proportional to jAm? |, respec- 
tively. Their relative amplitude(oscillation intensity), is 
about 40 : 2 : 1 from a global fit [l9| of mixing pa- 
rameters as listed in Table HI The oscillation component 
1 — P21 dominates the Pee oscillation, while P31 and P32, 
which are sensitive to the neutrino mass hierarchy, are 
suppressed by the small value of sin^(26'i3). 

The observed neutrino spectrum in L/E space, taking 
the baseline L to be 60 km and all the other parame- 
ters from Table U except sin^(20i3), is shown in FiglU 
together with that of no oscillation. For comparison, the 
oscillation spectrum without P31 and P32 are also shown. 
For a very small sin (2^13), a normal x^ analysis on the 
L/E spectrum with binned data, which requires accu- 
rate knowledge on the neutrino energy spectra and much 
smaller binning than the energy resolution, is difficult for 
the mass hierarchy study. 

Since neutrino masses all appear in the frequency do- 
main as shown in Eq. 21 a Fourier transform of F{L/E) 
shall enhance the sensitivity to the mass hierarchy. The 
frequency spectrum can be obtained by the following 
Fourier sine transform(FST) and Fourier cosine trans- 
form(FCT): 
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FIG. 1: Reactor neutrino spectra at a baseline of 60 km 
in L/E space for no oscillation (dashed dotted line), 1 — P21 
oscillation (dotted line) and P^e oscillation in the cases of NH 



and IH, assuming sin (26i 
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F{t) cos{Lut)dt 
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where uj is the frequency, uj — 2.54A?n^-; t 
variable in L/E space, varying from tmin 
t 



^ is the 
-J^ to 



max — E ' 

Since Pee is a linear combination of 1 — P21, P31 and 
P32, EST and FCT spectra can be divided into three 
components corresponding to 1 — P21, P31 and P32 re- 
spectively. FigO shows the three components of the EST 
and FCT spectra together with full Pee oscillation for 
both NH and IH cases. The oscillation frequency is pro- 
portional to Am? , so we can scale the frequency to be 
6m^ and plot the spectra in axis of Sm'^ in the interested 
frequency range of 1.8 xlO~^eV^ < Sm'^ < 3.0xlO~^eV^. 
From Eiglll we know that: 

1. P31 and P32 components dominate the FCT and 
FST spectra in the interested frequency range of 
l.SxlO-^eV^ < (5m2 < S.OxlO-^eV^ since |Ato|i| 
and |Ar7T,|2l are in this range, while 1 — P21 is very 
weak since its oscillation frequency is in a much 
lower range. The FST and FCT spectra of Pee are 
approximately the sum of P31 and P32 components 
which are sensitive to mass hierarchy. 



2. For NH, the P32 FCT and FST spectra are left- 
shifted with respect to the P31 spectra because 
< |Am|j|; while for IN, the P32 spectra 



\Aml^\ 



are right-shifted because |Ato§2| > |Am§;^|. 

3. The peak of FCT spectrum corresponds to the zero 
point of FST spectrum. This feature is helpful to 
identify the position of |Am§2 1 and | Ato|;^ |, without 
knowing their accurate values a priori. 

4. For FCT spectrum, P32 and P31 components have 
similar shapes with the peak around |Am|2| and 
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FIG. 2; Fourier sine (FST) and cosine (FCT) transform 
spectra for 1 — P21 component (dotted line), P32 component 
(dashed line), P31 component (dot-dashed line) and all the 
components of Pee (solid line) in the cases of NH and IH. 
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respectively, and two valleys on each side 
of the peak. The amplitude of P32 to that of P31 has 
a ratio of about 1:2 determined by tan^(6'i2). The 
shapes of P32 and P31 are left-right symmetric with 
respect to their peaks (mirror symmetric). This 
symmetry is broken for Pee as an approximate sum 
of P32 and P31 in different ways for NH and IH. For 
NH, the peak of P32 is at the left of the valley of 
P31, while for IH, the peak of P32 is at the right 
of the valley of P31. This feature can be used to 
distinguish NH and IH. 

5. For FST spectrum, the shapes of P32 and P31 are 
positive-negative symmetric with respect to zero 
(rotation symmetric) around |Ato|2| and |A77i§]^|, 
respectively. This symmetry is broken for Pee in 
different ways for NH and IH. For NH, the peak of 
P32 is at the valley position of P31 , while for IH, the 
valley of P32 is at the peak position of P31. This 
feature can be also used to distinguish NH and IH. 

As discussed above and shown in Fig. [2l the normal 
or inverted mass hierarchy can be distinguished by the 
symmetry breaking features of the FCT and FST spectra. 
To quantify these features, two parameters, RL and PV, 
are introduced as the following: 



RL 



RV -LV 
RV + LV 



, PV 



P-V 

P + V 



(6) 



where RV is the amplitude of the right valley and LV is 
the amplitude of the left valley in the FCT spectrum. P 



is the amplitude of the peak and V is the amplitude of the 
valley in the FST spectrum. From the above discussion, 
we know 



RL>0 a.nd PV >0 ^ NH 
RL<0 and PV <0 ^ IH 



(7) 



The values of RL and PV as well as the shapes of FCT 
and FST spectra depend on the baseline and neutrino 
mixing parameters. Parameters such as sin 612, Arn^i, 
and Arn^2 are relatively well known, hence only small 
uncertainties are introduced. The baseline and sin (2^13) 
are more important and are discussed below. 

1. Baseline determines the oscillation cycles. To max- 
imize the symmetry breaking of FCT and FST 
spectra, we scan the baseline length and find that 
the peak (valley)of P32 spectrum lays on the val- 
ley (peak) of P31 spectrum around 60 km. The 
widthes of peaks and valleys of the Fourier spectra, 
which are proportional to 1/L, are also determined 
by baseline. In an extreme case, the peak and val- 
ley of P31 and P32 spectra all become (5-functions at 
infinite baseline, hence are well separated from each 
other. In fact, this is already the case at 200 km 
and the mass hierarchy can be determined by look- 
ing at the position of the smaller peak (P32 compo- 
nent). If it is on the left side of the main peak (P31 
component), it is NH. Otherwise it is IH. However, 
since the neutrino flux from reactors is proportional 
to l/L"^, shorter baseline, say at 60 km, is the best 
from an experimental point of view. The actual op- 
timum baseline can be determined by taking into 
account both statistical and systematical errors. 

2. sin^(20i3) determines the amplitude of the Fourier 
spectra of P31 and P32. At sin^ (26*13) = 0, P31 
and P32 components will vanish and no features 
can be used to discriminate the mass hierarchy. A 
minimum value of sin^ (26*13) to distinguish NH and 
IH experimentally will be analyzed by taking into 
account possible experimental errors [20[. 

In order to understand the robustness of the discrim- 
ination method using FCT and FST spectra, values of 
baseline are scanned from 46 to 72 km; sin^(20i3) from 
0.005 to 0.05. The resultant RL and PV values are weU 
separated into two clusters, corresponding to the case of 
NH and IH respectively, as shown in FigO 

The FCT and FST spectra for sin2(26li3) = 0.005 are 
shown in FigH) Although a detailed experimental anal- 
ysis of error contour is to be completed [20] , the features 
of NH and IH are still very distinctive. On the FCT spec- 
trum, a valley appears at the left of the prominent peak 
for IH, and a peak appears at the left of the valley for NH. 
On the FST spectrum, there is a clear valley for IH, while 
for NH it is a peak. In comparison, the Fourier power 
spectrum used in Ref. [l^ is also shown in Figdl The 
FCT and FST method is more sensitive than the Fourier 
power spectrum method for a very small sin (2^13). 
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FIG. 3: Distribution of RL and PV values for different pa- 
rameters of baseline and sin^(2Si3). For each parameter to 
be scanned, the default baseline is 60 km and all the other 
parameters are the values as in Table HI Two clusters of RL 
and PV values are clearly seen for NH and IH cases. 



For even smaller sin^(26'i3), the main peak becomes 
less significant. For example, if the main peak is required 
to be twice higher than that of noise, siv?{29i^) must be 
greater than 0.005 in order to clearly identify the main 
peak, for a variety of neutrino energy spectra in a rea- 
sonable range. 

For a reaUstic experiment in the near future, the en- 
ergy resolution and statistics are of the most concern. At 
60 km, 9i2 has the least impact to the mass hierarchy de- 
termination. The energy resolution must be good enough 
not to smear the difference between P31 and P32, which 
requires the energy resolution be better than 3%/-\/!e. A 
detector with a mass at 10 kton level may be necessary, 
depending on the size of ^13. If shortening the baseline, 
the noise in the Fourier spectra from 612 oscillation in- 
creases, thus degrade the sensitivity. In the mean time 
requirements to the energy resolution and the detector 
size are relaxed. The optimization of the baseline as well 
as the energy resolution and detector size for different 6113 
assumptions are undergoing. 

In summary, the method to discriminate the mass hier- 
archy has been studied by using a Fourier sine(FST) and 
cosine (FCT) transform to the observed reactor neutrino 
L/E spectra. The FCT and FST spectra can separate P31 
and P32 oscillation components from the large 1 — P21 
component in a specific 5m? range. Features of mass 
hierarchy are enhanced in this representation and more 
sensitive than that of the Fourier power spectrum at very 
small sin^(26'i3). Wc found that an ideal detector at an 
intermediate baseline (~ 60 km) could identify the mass 



hierarchy for a mixing angle sin (26*13) > 0.005, without 
requirements on accurate information of reactor neutrino 
spectra and the value of A?ti|2- A paper of a detailed 
analysis of experimental errors will be released soon [20| . 
Similar methods can be applied to other experiments us- 
ing different neutrino sources, such as accelerator-based 
neutrino beams or atmospheric neutrinos. 
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FIG. 4: The FCT and FST spectra and Fourier power spec- 
trum for sin^(26li3) = 0.005. The sohd hue is for NH and the 
dashed line is for IH. The FCT and FST spectra have distinc- 
tive features to identify the mass hierarchy, which looks more 
sensitive than the Fourier power spectrum method. 
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